Diffusion on a stepped substrate 



Juha Merikoskiu and See-Chen Ying 
Department of Physics, Brown University, Box 1843, Providence R.I. 02912 

(February 1, 2008) 



0\ 



O 



> 
00 

o 



o 



We present results for collective diffusion of adatoms on a stepped substrate with a submonolayer 
coverage. We study the combined effect of the additional binding at step edge, the Schwoebel barrier, 
the enhanced diffusion along step edges, and the finite coverage on diffusion as a function of step 
density. In particular, we examine the crossover from step-dominated diffusion at high step density 
to terrace-dominated behavior at low step density in a lattice-gas model using analytical Green's 
function techniques and Monte Carlo simulations. The influence of steps on diffusion is shown to 
be more pronounced than previously anticipated. 

PACS numbers; 68.35.Fx, 05.40. +j 



The importance of surface steps in dynamical processes 
such as growth and ordering is widely recognized [Q. 
Among other effects, there exist extra activation barriers 
at the steps leading to a different mobility for adsorbates 
there as compared to the flat terrace region In re- 
cent years, several new techniques have been employed to 
measure the diffusion rate of adsorbate on well-defined 
low index surfaces |^-^. However, since most of these 
experiments involve studies of density fluctuations over 
a length scale much larger than the average distance be- 
tween steps on the surface, it is important to ascertain 
whether the measured diffusion constant corresponds to 
the intrinsic value for the wide terrace region, or whether 
it is strongly influenced by the rates near the step edges. 
While many theoretical studies on the extra barriers at 
step edges exist, only a few consider the effect of steps 
on measurable diffusion constants and criteria for 

the crossover from step-dominated to terrace-dominated 
diffusion with increasing terrace width are not well estab- 
lished 0. 

In this Letter, we present a study of the effect of vari- 
ous microscopic modifications of adatom jump rates due 
to surface steps on macroscopically observable collective 
diffusion. We describe the system of adsorbates on a 
substrate with a periodic array of steps by a lattice-gas 
model. We use a combination of analytical approach and 
Monte Carlo simulations to analyze the dependence of 
diffusion on step density in this model. We discuss the 
implications of our results on the experimental determi- 
nation of diffusion constants, and the recent measure- 
ments by Xiao et al. [^-|6| using optical diffraction tech- 
niques in particular. 

Following an earlier work dealing with the zero cov- 
erage limit 0, we introduce several energy barriers in 
our model characterizing the adsorption on the stepped 
substrate. These are the additional binding energy Eb 
at lower step edge, the Schwoebel barrier Eg for motion 
across step edges, the diffusion barrier on terraces -Eqj 



and the barrier for jumps along lower step edges E2 
The corresponding potential profile in direction perpen- 
dicular to step edges [x direction) is shown in Fig. 0(a). 
These energy barriers lead to the following jump rates as 
shown in Fig. |l|(b) and Fig. ^c) 



Fo = e-^"/^^ = Fie^-/'^^ 
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The substrate surface has a periodic array of steps sep- 
arated by terraces of width L lattice sites [see Fig. |^(c)]. 
Thus each adsorption site in the lattice gas model can be 
labelled by the coordinate (x, y) of the unit cell together 
with a site index s = 1,2, ... ,L within the unit cell. We 
assume hard-core interactions between the adatoms. For 
each lattice site we then define a stochastic occupation 
variable n%y{t), which can take only the values and 
f . To analyze collective diffusion, we define the Green's 
functions GlCy,{t) = -2TTie{t){u%y{t)u%, where 
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2, . . .L, the fluctuation variables 



s,s' = l, , 

fined by <y (t) = <y (t) - « 
step function ||9|-pl|. 

The equations of motion for G%y^,y 



%y{t) are de- 
), and 9{t) is the Heaviside 



[t) follow from the 



rate equations for n^^^ |12| . After Fourier transforming 
with respect to the time and cell coordinates, we can 
obtain a closed set of equations for the Green's func- 
tions (w) in (k, w) space by decoupHng the higher 
order Green's functions. In this mean-field approxima- 
tion, Gl^ {lo) turn out to be inversely proportional to 
the determinant of an almost tridiagonal L x L matrix 
yl(k, a;). The interior elements of A, 



Ass = iiu + 2To cos kyb — ATq 

As,s±l — Fq, 



(2) 



are related to terrace jumps, and the modified elements 
at the corners of A, 
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All = iLo+2T2coskyh -Ti-Tu-2T2 -(ro-ri+rd-r„)c( 

A22 ^iuj + 2Tq cos kyh - 4ro - {Ti - ro)ce 
All =iuj + 2ro cos kyh - SFo - - (r„ - Td)ce 
A12 = To + (Fi - ro)ce 
A21 - Ti + (Fo - ri)ct 
Ail = [r<i + (r„ - rrf)ce]e-'^'=^'^ 

Ail = [r„ + {Yd - ^„)Q]e+*^'=^^ (3) 

are related to jumps near step edges. All the other el- 
ements of A are identically zero. Here a and 6 are the 
nearest-neighbour distances in x and y directions, respec- 
tively. By symmetry, the average occupation numbers 
{n%y) depend only on the row index s. In the present 
model, we have only two distinct row coverages Ce and 
Cf, defined by Ce — {n^y) for the lower edge row and 
Ct = (n'^y) = . . . = {n^y) for terraces. The value of total 
coverage c — ct + {Ce — Ct)/L together with the detailed 
balance condition 



Ce(l - Ct) 
Ct(l - Ce) 



(4) 



completely determine the values of Cg and ct- 

The diffusion tensor D can be extracted from the dif- 
fusive poles of the Green's functions, i.e. from the zeros 
of det A in the hydrodynamic limit k and lu 
at uj = — ik-D-k. By expanding the determinant and 
collecting the leading terms for small uj and k, we find 
the diagonal elements of D to be [T^ 



(l-c,)(l-ct)ctroLV 



A, 



K [{L - 2)(1 - Ce) + (1 - Ct){l + eEs/kT)^ 
-[{L - 1)(1 - ct)ctTo + (1 - Ce)cer2]b^ 



K= (i-l)(l-Q)ct + (l-Ce)Ce. 



(5) 



By symmetry, the nondiagonal elements of D are identi- 
cally zero. 

To test the mean-field result of Eq. ||, we have per- 
formed standard Monte Carlo (MC) simulations for 
the same model. In the MC simulations, the collective 
diffusion constant has been calculated using the Green- 
Kubo relation [01 



D 



uv — f lim 1 



N 



^[r«(i)-r«(0)] 
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where the sum is over all particles and v = x,y. For our 
model, the 'thermodynamical factor' /, which is inversely 
proportional to particle number fluctuations jsj, can be 
calculated exactly yielding the result / = Lc/ k with k 
deflned in Eq. (jS) The coverage dependence of D 

as given by Eq. (m together with the simulation results 
are shown in Fig.n2 for a temperature comparable to the 
various activation barriers in the system: kT = 1/2 with 



Eb = Es = Eq — E2 = 1- Typical simulation param- 
eters were 10® MC steps per atom for thousand atoms 
in the system. We conclude that the mean-field theory 
agrees very well with the MC result, with no discernible 
systematic deviation |l^. At low temperatures / has 
a maximum at c « 1/L due to the suppression of fluc- 
tuations when lower edge rows get filled (cg « 1) and 
terraces are almost empty (ct « 0), resulting in a maxi- 
mum of Dxx at c « 1/L. In the case of Dyy, this effect 
of reduced fluctuations is cancelled out by the enhanced 
blocking of jumps along the lower step edges. 

We now focus on the dependence of the diffusion con- 
stant on the width of the terrace as given by the expres- 
sion in Eq. (||) . In this regard, Dyy has a relatively weak 
dependence on the width, since diffusion in the direction 
parallel to the step edge has to proceed through the ter- 
race region as soon as the step edge row is filled. Below 
we shall concentrate on D^x, i-e. diffusion perpendicular 
to the steps. At very high temperature or for very wide 
terraces, the diffusion constant D^x would approach the 
value Doo = Yoa^ = e^^^^^^'^ a? appropriate for a single 
terrace, and the effect of the steps vanishes. In the other 
limit of low temperature and narrow terrace width, dif- 
fusion is dominated by the steps, and D^x oc e^'"^/^'^ ^ 
where the total barrier Esd is now Eq + Es + Eb ■ The 
crucial question is how the system crosses over from one 
limit to another. The crossover behaviour can be under- 
stood most easily in the limit e.^slkT ^ ^ and c ^ 1/i 
p^ , in which the expression of D^x in Eq. (||) simplifies 
to the form 



i + c(l + e-Bs/fcT)g£;B/fcT^ 



(7) 



from which we see immediately that the crossover occurs 
at L « ce^EB+Es)lkT ^ Pqj. ^^^^^ vtAvsis of {EB + Es)/kT, 

this implies that the effect of steps could be substan- 
tial even for terrace widths up to several hundred lattice 
spacings. 

Recently the diffusion of CO on clean Ni(llO) was stud- 
ied by Xiao et al. [^j^ . From a measurement on a sample 
with high step densities [|| , they deduce the value of the 
step-dominated barrier Esd for diffusion to be 239 meV. 
The activation barrier for CO diffusion along (001) on a 
good sample with average step separation of L w 170 was 
determined to be 120 meV. For diffusion along (110), a 
barrier of 95 meV was measured. These authors then use 
a heuristic argument to arrive at the following crossover 
criterion: « e'^EB+Es)/kT ^ Based on this formula they 
conclude that in the experimental temperature range of 
100 to 200 K, the influence of the steps should be negligi- 
ble [|. The measured value of 120 meV is then assigned 
to be the value of the barrier Eq appropriate for CO dif- 
fusion along (001) on clean Ni(llO). 

We now apply our results to analyze whether the effect 
of steps is indeed negligible under these conditions in this 
system. We take as input the measured value for iJsd of 
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239 meV and the estimated values oi L = 170, Eq = 120 
meV for the (001) direction and L = 520, Eq — 95 meV 
for the (110) direction, and set Es — |]l6|. For simphc- 
ity, we have also taken all the prefactors of the various 
jump rates to be equal |l^. We calculate the effective 
extra barrier due to the steps, AE{T,L), for diffusion 
perpendicular to the steps, defined as the local slope of 
the Arrhenius plot of D^x/Doo using Eq. (|) for cover- 
age c — 0.5. In Fig. H we show the result as a function 
of temperature for diffusion along the two directions. In 
both cases, below 300 K, AE is sizable and the effective 
barrier is not just the terrace value Eq but considerably 
affected by steps. Clearly the heuristic crossover criterion 
in Ref. underestimates the influence of the steps. 

In a more recent experiment |^, the effect of sulphur 
as an impurity on diffusion of CO on Ni(llO) was stud- 
ied. It was proposed that the main effect of sulphur is to 
increase the effective step barrier dramatically. A barrier 
of 323 meV was measured for diffusion along both the 
main orientations. These authors conclude that in this 
case, the diffusion is dominated by the sulphur-poisoned 
step edges, and the measured value of 323 meV is then 
assigned to be £^sd, i-e. the step-dominated barrier. We 
have also applied our analysis to this system and find 
that here with L = 520 assuming Es = 0, the steps in- 
deed dominate for Eq < 200 meV. So the determination 
of the sulphur-modified step barrier in Ref. is rea- 
sonably consistent with the starting assumption in this 
case. 

In this Letter, we have presented a theoretical analysis 
of collective diffusion on stepped substrates, and studied 
the dependence of the diffusion constant on the terrace 
width, coverage and temperature. In particular, we have 
established a criterion for crossover from terrace domi- 
nated diffusion to step dominated diffusion. Compared 
with a heuristic criterion used before our theory 

predicts a much stronger influence of the steps on the 
diffusion barrier. The experiments of CO on Ni(llO) in 
Refs. [^-|| represent a situation characteristic of diffusion 
on smooth metal surfaces with steps as far as few hun- 
dred lattice spacings apart and temperature range of 100 
to 200 K. Our results show that under these conditions, 
the influence of the steps on the measured diffusion bar- 
rier can be considerable, and a careful theoretical analy- 
sis is needed to interprete the measured barrier value. In 
this study, the only interaction effect taken into account 
is the exclusion of double occupation on the same site. 
The effect of additional interaction between adatoms on 
the influence of the steps remains to be investigated. 
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FIG. 1. Geometry and jump rates of the lattice model for 
diffusion on a stepped substrate, (a) The potential profile in 
X direction, (b) Side view of the model showing the various 
hopping rates for jumps in x direction near step edges, (c) 
Top view of the model showing the size of one unit cell with 
the indices s = 1, 2, . . . , L of each lattice site within the cell. 
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FIG. 2. Mean field and Monte Carlo results for Dxx and 
Dyy as a function of coverage at temperature kT = 1/2 for 
the hard-core model with Es = Eb — Eq — E2 = 1. Lines 
show the mean field solution of Eq. (jsj) , and plotting symbols 
denote the results of Monte Carlo simulations with errorbars 
less than the size of the symbols. Here diffusion coefficients 
are shown in units of Doo that is the value for infinite terrace 
width and coverage in units of one monolayer. 
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FIG. 3. The effective extra activation barrier due to the 
steps for diffusion of CO along the (001) and (110) directions 
on Ni(llO). See text for parametrization. 
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